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Multiple-Retrapping Process in High-T c Intrinsic Josephson Junctions 
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We report measurements of switching-current distribution (SWCD) from a phase-diffusion branch 
to a quasiparticle-tunneling branch as a function of temperature in a cuprate-based intrinsic Joseph- 
son junction. Contrary to the thermal-activation model, the width of the SWCD increases with 
decreasing temperature, down to 1.5 K. Based on the multiple-retrapping model, we quantitatively 
demonstrate that the quality factor of the junction in the phase-diffusion regime determines the 
observed temperature dependence of the SWCD. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.40.+k, 74.45. +c 
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The escape of a system trapped in a metastable state 
governs the reaction rate in various dynamical systems, 
where the escape is made by a noise-assisted hopping [l| . 
In the case of a Josephson junction (JJ), a thermal noise 
induces an escape of a phase particle representing the sys- 
tem from a local minimum of the potential well. In an 
underdamped JJ with hysteresis in the voltage-current 
( V-I) characteristics a single escaping event induces a 
switching from a zero-voltage phase-trapped state to a 
finite- voltage phase-rolling state [2J. In an overdamped 
JJ without hysteresis, however, the energy of an escaped 
phase particle is dissipated during its motion so that the 
particle is retrapped in another local minimum of the 
potential. The phase particle repeats this thermally ac- 
tivated escape and retrapping process, i.e., the multiple- 
retrapping process. It results in a phase-diffusion branch 
(PDB); a resistive branch with a small but finite voltage 
for a bias current below a switching current, Isw [3||- 

A hysteretic JJ can also evolve into this multiple- 
retrapping state with increasing temperature (T) when 
the energy fed by a bias current, near a mean switch- 
ing current that is suppressed by an increment in T, 
gets comparable to the dissipated energy. Recently, this 
phase-retrapping phenomenon in the hysteretic JJ has 
been intensively studied in association with the T depen- 
dence of a switching current distribution (SWCD) [4], [5(. 
The main finding of these studies is that the retrapping 
process in the hysteretic JJ modifies the switching dy- 
namics, in such a way as to reduce the width of SWCD 
with increasing T. This SWCD behavior, contrasting to 
the usual thermal activation model, has also been sug- 
gested to be caused by the enhanced dissipation due to 
a high-frequency impedance, which originates from the 
contribution of the impedance of the measurement lines 
[a, ll| , incorporating into temperature effect [8( . At suffi- 
ciently high T, the hysteretic JJ begins to show a PDB. 
This PDB in a hysteretic JJ can be found whenever the 
thermal energy, fcgT, is comparable to the Josephson en- 
ergy, Ej (=hl c /2e; I c is the noise-free critical current 
and ks is the Boltzmann constant). Thus, an ultra- 
small hysteretic JJ exhibits the PDB even at a-few-mK 



region [9(. As the dissipation is represented by the qual- 
ity factor Q of the JJ, however, it has not been clearly 
resolved yet how Q, depending both on the temperature 
and the impedance, affects switching between a PDB and 
a quasiparticle-tunneling branch (QTB). 

In this letter, we study switching from the PDB 
to the QTB in intrinsic Josephson junctions (IJJs) of 
Bi2Sr2CaCu208+ K (Bi-2212) at various temperatures. 
The effective junction impedance Zj including the in- 
fluence of an external circuit in a PDB, was estimated by 
fitting V-I curves to the phase-diffusion model. The val- 
ues of Zj turn out to be an order of the measurement-line 
impedance Zl, which indicates that the phase dynam- 
ics of our IJJs is governed by the environmental dissi- 
pation. Our measurements of the T dependence of the 
switching rate Ts and the corresponding SWCD in a 
single IJJ are in good agreement with those estimated 
by the multiple-retrapping model. This study quanti- 
tatively clarifies how the impedance- and temperature- 
related dissipations, and the corresponding Q determine 
the switching dynamics in a hysteretic IJJ with phase- 
diffusion characteristics. 

A stack with the lateral size of 2.5x2.9 /xm 2 in Bi- 
2212 single crystal was defined by using focused-ion-beam 
(FIB) process [10j. High- intensity FIB irradiation is 
known to degrade theperipheral region by the scattered 
secondary ion beam [ill ]. In the milling process, we used 
a relatively high ion-beam current of 3 nA, correspond- 
ing to ~200 pA//im 2 . This high ion-beam current much 
reduced the interlayer tunneling critical current density 
down to ~8 A/cm 2 in A^=12 junctions out of total ^100 
junctions in a stack, which were estimated from the num- 
ber of QTBs in V-I curves of the stack. Four-terminal 
transport measurements were carried out in a pumped 
He 4 dewar with the base temperature of 1.45 K. Room- 
temperature 7r-filters were employed and measurement 
lines were embedded in silver paste at cryogenic tem- 
peratures to suppress high-frequency noises propagating 
along the leads. The measurements were made by using 
battery-operated low-noise amplifiers (PAR-113). The 
ramping speed of the bias current and the threshold volt- 
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FIG. 1: (color online) V-I characteristics at various temper- 
atures of the Bi-2212 IJJs stack. The upward and downward 
arrows indicate the switching and the return currents, respec- 
tively, at T=15.1 K. Black solid curves are the best fits to the 
phase diffusion model at three different temperatures of 15.1 
K, 18 K, and 23 K, with the effective junction impedance 
Zj — 125, 125, and 128 Q as the best-fit parameters, respec- 
tively. Inset: the R-vs-T curve, where T c is indicated by an 
arrow. 



age (obtained from the maximum voltage of the PDB) in 
measuring I sw were 7=30 ^iA/sec and 14^=110 /xV [12j |. 
respectively. Measurements were made on the first jump, 
which connected the PDB to the first QTB of the weak- 
est junction in the stack. For each distribution, 10000 
switching events were recorded with the current resolu- 
tion of 90 pA. 

The inset of Fig. 1 shows the tunneling R vs T curves 
with the superconducting transition at T c =90 K, which is 
indicated by an arrow. The c-axis tunneling nature in the 
Bi-2212 stacked junctions is evident by the increasing re- 
sistance with decreasing T above T c . At T<T C , the junc- 
tion resistance does not vanish completely, even down to 
T=10 K, which is attributed to the phase diffusion in the 
junctions. Fig. 1 shows the V-I curves at various T be- 
low T c . The nonlinear curve at T=50 K begins to show 
a hysteresis below T~30 K. For a bias below a Isw in 
each curve, a pronounced non-zero resistive branch ap- 
pears below T=23 K. Here, the finite resistance is caused 
by the phase diffusion in the condition of ksT~2Ej [3J. 

J predicts the current- 
T ^ with 



The phase-diffusion model [13|, ll 
voltage relation of I(V) = AeIswk 



h v , -> + [2eZjk B T/h] 2 

a frequency-dependent effective junction impedance Zj. 
Here, the phase diffusion is assumed to take place in all 
twelve junctions with FIB-suppressed critical currents, so 
that the voltage bias per junction becomes v = V/N . We 
use this relation to fit the PDB as shown by black curves 
in Fig. 1 with Zj as the best-fit parameter. Zj turns 
out to be 125, 125, and 128 fi for T=15.1, 18, and 23 K, 
respectively. These values of Zj are comparable to the 
measurement line impedance Z^ (50^100 Q,), even with a 
much higher quasiparticle tunneling resistance of i?„ p ~35 



kfi as seen in Figure 1. This implies that the dissipation 
in the PDB is dominated by the high-frequency (of an 
order of the Josephson plasma frequency, u p ) dissipation 
through the measurement lines [y, |2fl . Since the energy 
of the escaped phase particle is dissipated through the 
environment in the phase-diffusive regime, the switching 
also becomes sensitive to this dissipation process. 

Now, we turn to the switching event from the PDB to 
the QTB. Fig. 2(a) shows the SWCD (scattered sym- 
bols) at various temperatures. The standard deviation 
(a) and the mean switching currents ({Isw)) are shown 
as a function of T in the inset of Fig. 2(a). The T 
dependence of a contradicts to that of a conventional 
underdamped J J, where a increases with T in a thermal- 
activation regime [2j. Fig. 2(b) shows Ts (scattered 
points) vs I calculated from the SWCD of Fig. 2(a) fol- 
lowing the Fulton and Dunkleberger analysis [2j. With 
lowering T, Ts{I) shows a pronounced change from an 
almost linear to a down-turn nonlinear bias-current de- 
pendence in a semi- logarithmic plot. To explain this be- 
havior, we adopted the multiple-retrapping model [5j, |gfl . 
In the phase-diffusive regime, the successive retrapping 
processes suppress the switching rate, Ts- The switching 
to the high-voltage branch (i.e. the QTB) occurs only 
when the phase particle is not retrapped after escaping 
from a local energy minimum. A phase particle escaped 
from a potential minimum has a probability, Prt, to be 
retrapped in the next potential minimum. The switching 
rate Ts, including Prt, is expressed by [8(] 



Ts — Tta(1 — Prt) 



ln(l - Pj 



RT) 



P 



RT 



(1) 



Here, r Tj 4[=^exp(--A^L)] is the thermally acti- 
vated escape rate, w p =w p o(l-7 2 ) 1 / 4 , u} p o={2eI c /hC) 1 / 2 , 
A[/(</>)[=2.Ej((l-7 2 ) 1 / 2 -7arccos7)] is the escape energy 
barrier, and 7[=J/i c ] is a normalized bias current. 
The retrapping probability can be obtained by an 
integration (see Ref. |15| ) of the retrapping rate 
r ^-27c(AC/ fl T/7rfc B r) 1 / 2 exp(-AC/ flT /A :B T), where 
/^URT(I)—Z 2 jC(I-I r o) 2 and I r o is the noise-free return 
current from a QTB to a PDB [161] , Fig. 3(a) shows the 
experimental switching distribution (red dots) at T=1.5 
K with the corresponding fit (green curve) obtained by 
using Eq. (1) with the best-fit parameters of Zj— 61.9 O, 
7 C =1.26 /iA, and J r o=63 nA. The junction capacitance, 
330 fF, was estimated from the typical value of 45 fF//im 2 
for Bi-2212 IJJs [12J. The corresponding switching rate 
(red dots) and the fit (green curve) are shown in Fig. 
3(b), with the same parameters. An excellent agreement 
is obtained in both fittings. 

These results are analyzed in terms of Zj- and In- 
dependence of Prt- The inset of Fig. 3(a) shows the 
calculated Prt-vs-I curves for various Zj at T=1.5 K, 
with the values of I c and I r o obtained from the best fits 
of Figure 3. The retrapping-probability curve shifts to 
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FIG. 2: (color online) (a) The switching current distribution 
and (b) the corresponding Ts at various T from 1.5 K to 4.2 K. 
The solid lines are calculated based on the multiple-retrapping 
model. Inset: Temperature dependence of the standard de- 
viation (a) and the mean switching current ({Isw))- Ar is 
the range of the observable switching rate for a given If,. In- 
set: the retrapping current distribution (solid circles) and the 
calculated behavior (solid line) at T=4.2 K. 



higher currents as Zj decreases, due to the presence of 
a Zj dependence of AUrt in the exponential factor of 
Tut- The current positions of almost vanishing Prt, in- 
dicated by downward arrows, approximately correspond 
to the maximum current allowing the retrapping. We 
denote this current as Ipd- Physically this is the same 
current as the one denoted I m in Ref. [6j. The system 
can hardly be retrapped at a current higher than Ipo, 
because in this case the energy fed to the system by the 
bias current gets larger than the dissipated energy. By 
equating the energy fed and the energy dissipated, sim- 
ilar to McCumber and Stewart analysis [4J, one obtains 
the relation 



I PD = 4:I c /ttQ pd , 



(2) 



where Qpd is the phase-diffusion quality factor at uj^uj p 
[6J]. In fact, the noise-free retrapping current can be 
written in a form similar to Eq. (2), namely as I r Q = 
4I c /irQ(tu=0) 3]. The thick black curves in Figs. 3(a) 
and 3(b) show SWCD and the corresponding Tta(I), re- 
spectively Other solid curves in Fig. 3(b) are Ts(I) in 
Eq. (1) for varying Zj under the multiple retrapping pro- 
cesses. Ts{I) with each Zj in Fig. 3(b) starts to drop 
quickly from Tta(I) at 1=1 pd, which are denoted by 
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FIG. 3: (color online) (a) The solid lines represent the cal- 
culated switching current distributions for various Zj values 
of 106, 61.9, 50, and 40 Q, corresponding to Q PO =3.88, 2.37, 
1.90 and 1.55, respectively. Red dots show the SWCD at 
T—1.5 K. Inset: Theoretical retrapping probabilities versus 
bias current, for various Zj values at T—1.5 K. The down- 
ward arrows indicate the Ipd values, (b) Estimated switch- 
ing rate versus the bias current. Red dots correspond to the 
experimental switching rate for T—1.5 K. Estimated Ts for 
various Zj values in Eq. (1) are shown as solid lines. Inset: 
Estimated retrapping probabilities versus bias current with 
Zj=61.9 Q at various T. Thick black curves in (a) and (b) 
show the thermally activated switching current distribution 
and the corresponding escape rate, respectively, without re- 
trapping. 



arrows in the figure. We define Ipp, as the current value 
corresponding to Prt=0.01, where Vs(Jpd) is nearly the 
same as the Tta(Ipd) as shown in Fig. 3(b) [8(. The 
impedance of 106 £1 gives the same SWCD as thermally 
activated SWCD without retrapping because Ts(I) in 
Ar well overlaps with Tta(I) although the impedance 
is of an order of Zl- When T$ {Ipd) crosses over the 
bottom line of Ar while Zj keeps decreasing, the high- 
frequency dissipation affects the SCDW: the observable 
window of Ar (between the two dashed lines) for a fixed 
lb shifts to the steeper section with decreasing Zj, result- 
ing in the decreasing width of the SWCD in Figure 3(a) 
at a constant T. This shift itself is caused by a reduction 
of Qpd with decreasing Zj at a constant T. 

Since Ipp, is also affected by a change in T as fol- 
lowings, the shape of Ts(I) dose not simply follow the 
variation of Zj as T changes. The inset of Fig. 3(b) 



TABLE I: Fitting parameters for the switching events for se- 
lected temperatures 



T(K) 


/c(mA) 


J r0 (nA) 


Zj(n) 


I PD (nA) 


Qpd 


1.5 


1.263 


63.00 


61.9 


678 


2.37 


2.4 


1.209 


65.10 


77.6 


683 


2.25 


3.0 


1.006 


63.44 


86.0 


685 


1.87 


3.6 


0.735 


65.22 


94.7 


682 


1.37 


4.2 


0.572 


63.79 


101.5 


572 


1.27 



illustrates Prt vs I at various temperatures. Here, Zj 
is fixed at 61.9 Q and other parameters, except for T, 
are set to be the same as for the inset of Fig. 3(a). The 
current position of the zero-temperature curve, indicated 
by an upward arrow in the inset of Fig. 3(b), corre- 
sponds to the fluctuation- free return current I r o- The 
value of Ipd shown by arrows increases with increasing 
T. Thus, in effect, the retrapping probability curve shifts 
to higher currents as T is raised, due to the presence of 
a Boltzmann-type exponential factor in the expression of 
T RT . Fig. 2 illustrates the calculated SWCD and T S (I) 
as solid curves at various temperatures with the best-fit 
parameters listed in Table 1, which agree well with the 
observation. Here, since the dissipation effect by temper- 
ature is already determined by a bath temperature, the 
main fit-parameter becomes Zj in Figure 3(a). As shown 
in Fig. 2(b), with increasing T, the calculated Ts{I) in 
the window of AT shows steeper regions. It causes the 
reduction of the width of SWCD with decreasing tem- 
perature. The ratio of I c and Ipd becomes smaller with 
increasing T as in Table I. Ipo even equals to I c at T=4.2 
K. This behavior leads to the conclusion that Qpd de- 
creases with increasing T following Eq. (2) despite the 
increase of Zj with T [l8[. The slop of the calculated 
r<j(J) in the window of Ar at sufficiently high temper- 
atures is insensitive to T variations, which leads to the 
saturation of a with increasing T as the inset of Figure 
2(a). 

Now we consider the retrapping dynamics related to 
the transition from a QTB to a PDB in an underdamped 
J J. It can be shown that only the zero- frequency dissipa- 
tion plays an significant role in the retrapping dynamics 
[6j. The inset of Fig. 2(b) shows the stochastic nature of 
the retrapping current as shown by solid circles at T=4.2 
K. The calculated retrapping distribution (solid line) by 
the retrapping rate, Tpx, is consistent with the experi- 
mental result. The best-fit parameters are i r o=63.8 nA 
and Zj=10 kfl (^R qp ), where the estimated noise-free 
I t q matches with the value used in T$ fitting in Table 
I. The junction impedance Zj estimated from the return 
currents is significantly larger than Zj (lu^lu p ) for the ob- 
served switching events. It indicates that the retrapping 
phenomena from QTB to PDB are mainly determined by 
a zero-frequency damping with Q(0)=11.4 with / c =572 
nA and / r0 =63.8 nA at T=4.2 K as shown in Table I. 



In summary, we clearly show that the multiple retrap- 
ping processes in underdamped IJJs, with much sup- 
pressed critical current and high tunneling resistance, 
govern the switching from the PDB to the QTB. The pre- 
dicted SWCD and T$ in the multiple-retrapping model 
are in good agreement with the observed broadening of 
the distribution of switching currents with decreasing 
temperature. We also demonstrate that the change of 
the shapes of the observed SWCD and the Tg in vari- 
ous temperatures can be understood by impedance and 
temperature dependence of Qpd, taking the retrapping 
probability into account. As the macroscopic quantum 
tunneling has been observed recently in IJJs of Bi-2212 
single crystals [19J , this study provides in a quantitative 
manner the role of the dissipation in quantum devices 
based on cuprate-based JJs. 
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